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Abstract 

Over the last decade, combinatorial algorithms have been obtained for exactly solving 
several nonlinear convex programs. We first provide a formal context to this activity by 
introducing the notion of rational convex programs - this also enables us to identify a 
number of questions for further study. So far, such algorithms were obtained for total 
problems only. Our main contribution is developing the methodology for handling non- 
total problems, i.e., their associated convex programs may be infeasible for certain settings 
of the parameters. 

The specific problem we study pertains to a Nash bargaining game, called ADNB, 
which is derived from the linear case of the Arrow-Debreu market model. We reduce this 
game to computing an equilibrium in a new market model called flexible budget market, 
and we obtain primal-dual algorithms for determining feasibility, as well as giving a proof 
of infeasibility and finding an equilibrium. We give an application of our combinatorial 
algorithm for ADNB to an important "fair" throughput allocation problem on a wireless 
channel. 
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1 Introduction 



The fascinating question of computability of market equilibria, which has been studied exten- 
sively within TCS over the last decade, has provided the area of algorithms a new direction, 
namely the design of efficient combinatorial algorithms 1 for exactly solving nonlinear convex 
programs. For this purpose, the classical primal-dual paradigm was suitably extended from its 
usual setting of LP-duality theory to convex programs and the KKT conditions; this task was 
initiated in [DPSV08]. 

We note that all problems attacked so far were total, i.e., their convex programs always have 
finite optimal solutions. The main contribution of this paper is to develop the methodology 
for handling problems that are not guaranteed to always have a solution, i.e., their convex 
programs may be infeasible for certain settings of the parameters. It turns out that in the 
setting of LP's, the primal-dual algorithms for non-total problems are not more involved than 
those for total problems. We illustrate this by comparing the algorithms for the problems 
of maximum weight perfect matching and maximum weight matching in bipartite graphs in 
Section 10. However, the situation is quite different for convex programs; see Section 1.2 for a 
high level description and Section 10 for a detailed analysis. 

The specific question that led to these ideas was combinatorially solving Nash bargaining 
games. Nash bargaining [Nas50] is a central solution concept within game theory for "fair" 
allocation of utility among competing players in the presence of complete information; it has 
numerous applications and a large following, e.g., see [Kal85, TL89, OR94]. The solution to 
a Nash bargaining game is obtained by maximizing a concave function over a convex set, i.e., 
it is the solution to a convex program. If the conditions for efficiently running the ellipsoid 
algorithm hold [GLS88], in particular, if efficient separation oracles can be implemented for its 
constraints and objective function, then the solution can be obtained to any required degree of 
accuracy in polynomial time. Instead, we resorted to the design of combinatorial algorithms 
because they have several advantages over "continuous" algorithms, see Section 1.1. 

The specific Nash bargaining game we study in this paper is called the Arrow-Debreu Nash 
Bargaining Game, abbreviated ADNB. The setup is the same as the linear case of the Arrow- 
Debreu market model, but instead of resorting to the solution concept of a market equilibrium 
for reallocating goods among the agents, we resort to a Nash bargaining solution. Interestingly 
enough, our algorithm for solving ADNB reduces it to a new, natural market model, which 
we call the flexible budget market. 

In the next section, we first provide a formal, mathematical context to the new algorithmic 
activity mentioned above. This also enables us to identify a number of questions for further 
study, see Section 12. 

1 Informally, an algorithm that conducts a search over a discrete space. A key distinction is that whereas a 
continuous method gives a way of solving the LP or convex program underlying a given instance, and provides 
no insight about the problem itself, a combinatorial algorithm handles all instances of the specific problem by 
exploiting its special combinatorial structure, e.g., see [VazlOc] for a detailed argument. 
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1.1 Rational convex programs and combinatorial algorithms for them 

A combinatorial algorithm for exactly solving a convex program is possible only if it admits 
rational solutions. Starting with the classic Eisenberg-Gale convex program [EG59], whose 
solution yields an equilibrium for the linear case of Fisher's market model [DPSV08], in re- 
cent years, researchers have identified several interesting convex programs that always admit 
rational solutions. First, let us formally define this class. We will say that a nonlinear convex 
program is rational if, for any setting of its parameters to rational numbers such that it has a 
finite optimal solution, it admits an optimal solution that is rational and can be written using 
polynomially many bits in the number of bits needed to write all the parameters. 

This class turns out to be surprisingly rich - besides several Nash bargaining games [VazlOa], 
it captures the Fisher market model under several classes of utility functions: linear [EG59, 
DPSV08], utility functions defined via combinatorial problems, including some in Kelly's 
[Kel97] resource allocation model [JV08, CDV], spending constraint utilities [VazlOb, BDX10], 
and piecewise-linear concave utilities in a market model that allows for perfect price discrimi- 
nation [GV10]. It also captures the linear case of the Arrow-Debreu market model (this follows 
from the convex program of [Jai07] and a proof of rationality in [Eav76]). 

Of course, rationality does not guarantee the existence of a combinatorial polynomial time 
algorithm. However, it turns out that such algorithms are known for solving almost all of 
the rational convex programs mentioned above. The only exceptions are 2-agent markets and 
Nash bargaining games given in [CDV] and [VazlOa], and the linear case of Arrow-Debreu 
markets. For the former, polynomial time algorithms are given, in [CDV] and [VazlOa], using 
only an LP-solver. This leads to a tantalizing question: is the class of rational convex programs 
solvable in this manner, i.e., by polynomial time algorithms that use only an LP-solver? 

Both algorithmic approaches, continuous and combinatorial, are valuable in their own right 
and have different advantages. Indeed, it is the synergy between these two approaches that 
makes algorithm design a potent field. We point out the advantages of the latter approach 
below. 

We believe that the case for designing combinatorial algorithms for rational convex programs is 
as compelling as that for integral LP's, i.e., linear programs that always admit integral optimal 
solutions. The latter led to deep combinatorial insights and constitute a substantial part of 
the field of combinatorial optimization. In recent years, insights gained from combinatorial 
algorithms for convex programs have also led to major progress. For instance, the recent proof 
of membership in PPAD of markets under piecewise-linear concave utilities [VY10] followed 
from a new, combinatorial way of characterizing equilibria [DPSV08] and helped settle, together 
with [CDDT09, CT09], the long-standing open problem of determining the exact complexity 
of this key case of markets. 

Combinatorial algorithms also have several advantages over continuous algorithms in applica- 
tions. For instance, recently Nisan et. al., faced with the problem of designing an auction 
system for Google for TV ads, converged to a market equilibrium based method, after explor- 
ing several different options [NBC + 09]. As stated by Nisan [Nis09], the actual implementation 
of this algorithm was inspired by combinatorial market equilibrium algorithms, which in turn 
solve convex programs combinatorially. Indeed, the easy adaptability of combinatorial algo- 
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rithms to the special idiosyncrasies of an application often makes them the preferred method. 
In Section 4 we present an application of our combinatorial algorithm for ADNB to a through- 
put allocation problem on a wireless channel. 

Although our algorithm is not strongly polynomial, recent results provide a reason to believe 
that our methodology could lead to it. Using scaling techniques, Orlin [OrllO] extended the 
DPSV algorithm to a strongly polynomial one, and an obvious open question is to extend the 
rest of the algorithms for rational convex programs as well. 

1.2 Determining feasibility of a convex program 

A linear program is infeasible iff no point satisfies all its constraints. In the case of a convex 
program, infeasibility could arise because of a second reason, namely the objective function is 
undefined at each point that satisfies all constraints. The convex program for ADNB exhibits 
only the second type of infeasibility; it always has points satisfying all constraints. 

Even for the case of total problems, the primal-dual method operates in fundamentally different 
ways in the two settings of integral linear programs and rational convex programs; for details, 
see Section 4 in [DPSV08] or Section 10 of the present paper. This difference manifests itself 
even more acutely when we move to non-total problems, as shown in detail in Section 10. 

For solving ADNB, we give a special procedure, Stage I in Algorithm 1, that tests for fea- 
sibility. If the instance is found to be infeasible, this procedure finds a proof of infeasibility. 
Otherwise, it yields special prices, called feasible prices, which provide Stage II with the right 
starting point for finding a solution to the given instance. 

We give two different proofs of infeasibility. The first one uses an explicit dual of the convex 
program for ADNB, found by Devanur [DevlO]; for proving infeasibility, our procedure shows 
that this dual is unbounded. Often it is difficult to find an explicit dual of a convex program. 
For this reason, we give a second proof of infeasibility. We give an LP whose optimal solution 
is positive iff the given instance is feasible. Therefore, a dual feasible solution of value at most 
zero is a proof of infeasibility. Our procedure also finds such a proof of infeasibility. 

1.3 New insights into balanced flows 

Differences between complementary slackness conditions for an LP and the KKT conditions 
for a convex program (see Section 4 in [DPSV08]) leads to new difficulties in the latter setting. 
The new algorithmic idea of balanced flows, introduced in [DPSV08], helps overcome these 
difficulties. Although this notion yields the desired efficient algorithm, the use of the h norm 
in the potential function used for proving polynomial time termination makes the proofs quite 
difficult. [DPSV08] observe that the I2 norm can be dispensed with for defining balanced flows 
and ask, "Can a polynomial running time be established for the algorithm using the alternative 
definition, thereby dispensing with the I2 norm altogether? At present we see no way of doing 
this ... ." 

In Section 11 we answer this question in the negative by providing an infinite family of exam- 
ples in which the natural l\ norm-based potential function makes inverse exponentially small 
progress whereas the I2 norm-based potential function makes inverse polynomial progress. 
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We observe that the combinatorial algorithms of [VazlOb] and [GV10], which solve rational 
convex programs, also use the notion of balanced flows crucially. Balanced flows play an even 
more fundamental role in our work. As detailed in Section 8.4, they are used in several crucial 
ways in both stages of our algorithm. In addition, they are also used for defining the central 
notion of feasible prices. 

2 The Nash Bargaining Game 

An n-person Nash bargaining game consists of a pair (Af,c), where Af C R™ is a compact, 
convex set and c G Af. Set Af is the feasible set and its elements give utilities that the n players 
can simultaneously accrue. Point c is the disagreement point - it gives the utilities that the n 
players obtain if they decide not to cooperate. The set of n agents will be denoted by B and 
the agents will be numbered 1, 2, ... n. Game (A/", c) is said to be feasible if there is a point 
v £ Af such that \/i G B, Vi > Cj, and infeasible otherwise. 

The solution to a feasible game is the point v G Af that satisfies the following four axioms: 

1. Pareto optimality: No point in Af can weakly dominate v. 

2. Invariance under afRne transformations of utilities: If the utilities of any player 
are redefined by multiplying by a scalar and adding a constant, then the solution to the 
transformed game is obtained by applying these operations to the particular coordinate 
of v. 

3. Symmetry: If the players are renumbered, then it suffices to renumber the coordinates 
of v accordingly. 

4. Independence of irrelevant alternatives: If v is the solution for (Af, c), and <S C R" 

is a compact, convex set satisfying c G S and v G S C Af, then v is also the solution for 
(S,c). 

Via an elegant proof, Nash proved: 

Theorem 1 Nash [Nas50] If game (Af, c) is feasible then there is a unique point in Af satis- 
fying the axioms stated above. This is also the unique point that maximizes IIj g £(^ — cf), over 
all v G Af. 

Most papers in game theory assume that the given Nash bargaining game (Af, c) is feasible. 
However, in this paper, it will be more natural to not make this assumption and to determine 
this fact algorithmically. Thus, we can have one of 2 outcomes: 

Thus Nash's solution to his bargaining game involves maximizing a concave function over a 
convex domain, and is therefore the optimal solution to the following convex program. 

maximize log(vj — cf) (1) 

i€B 

subject to v G Af 
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As a consequence, if for a specific game, a separation oracle can be implemented in polynomial 
time, then using the ellipsoid algorithm one can get as good an approximation to the solution 
of this convex program as desired in time polynomial in the number of bits of accuracy needed 
[GLS88]. 

3 The Game ADNB 

The game ADNB, short for Arrow-Debreu Nash Bargaining game, which will be studied 
extensively, is derived from the linear case of the Arrow-Debreu model. We state the latter 
first. Let B = {1, 2, . . . , n} be a set of agents and G = {1, 2, . . . , g} be a set of divisible goods. 
We will assume w.l.o.g. that there is a unit amount of each good. Let Uij be the utility derived 
by agent i on receiving one unit of good j; w.l.o.g., we will assume that is integral. If x-ij 
is the amount of good j that agent i gets, for 1 < j < g, then she derives total utility 

v i(%) ^ ' UijX{j. 

Finally, we assume that each agent has an initial endowment of these goods; for each good, 
the total amount possessed by the agents is 1 unit. 

W.l.o.g. we may assume that each good is desired by at least one agent and each agent desires 
at least one good, i.e., 

Vj G G, 3i G B : > and \/i G B, 3j G G : > 0. 

If not, we can remove the good or the agent from consideration. 

The question is to find prices for these goods so that if each agent sells her entire initial endow- 
ment at these prices and uses the money to buy an optimal bundle of goods, the market clears 
exactly, i.e., there is no deficiency or surplus of any good. Such prices are called equilibrium 
prices. 

The Arrow-Debreu market model gives one mechanism by which the agents can redistribute 
goods to achieve higher utilities. Another mechanism is to view this setup as a Nash bargain- 
ing game as follows. For each i G B, let q denote the utility derived by agent i from her 
initial endowment. Regard this as agent i's disagreement utility and redistribute the goods in 
accordance with the Nash bargaining solution. 

We will define game ADNB in a slightly more general manner: we will assume that the 
disagreement utilities, Cj's, are arbitrary numbers specified in the particular instance. (As 
stated in the Introduction, we will not deal with the nonsymmetric extension of ADNB in 
this paper.) Clearly, the Nash bargaining solution is the optimal solution to the following 
convex program: 

maximize logfvj — (h) (2) 

subject to Vi G B : V{ = UijXij 
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VjGG: 5> -<l 
Vz GB, Vj € G : > 

The KKT conditions for this program are: 

(1) Vj G G : Pj > 0. 

(2) VjeG: Pj >0^ EieB^j 

(3) \/zGB, VjeG:p, >^. 

(4) Vi EB, Vj € G : x„ > =>Pj 

Theorem 2 Program (2) is a rational convex program. Moreover, if it is feasible, then the 
dual solution is unique. 

Proof : We will show that if for a setting of rational parameters, program (2) is feasible, 
then the x^'s and pj's are solutions to an LP and are therefore rational numbers that can be 
written using polynomially many bits. First, guess the Xjj's that are non-zero in the optimal 
solution to program (2). Because of the assumption made on the instance, each pj will be 
positive. 

The variables of the LP will be the non-zero x^-'s and for each good j, a new variable qj, 
which is supposed to represent 1/pj- The LP will have the following constraints: for each qj, 
there is one equation corresponding to the KKT condition (2), and for each nonzero Xij there 
is one equation corresponding to the KKT condition (4). In addition, the LP has inequality 
constraints corresponding to KKT condition (3), for each i G B and j G G. In all these 
constraints, Vi is replaced by J2jeG u ij x ij- Finally, it has non-negativity constraints for all 
Xjj's and qj's. It is easy to check that all constraints are linear. 

Since program (2) is feasible, so is the LP corresponding to the correct guess. The solution to 
this LP will satisfy all KKT conditions and hence is an optimal solution to program (2). This 
established rationality of the convex program. Finally, the strict concavity of the objective 
function of program (2) and the fact that the convex combination of any two of its feasible 
solutions is also feasible implies that the optimal values of v^s is unique. Now, using the KKT 
condition (4), we get the uniqueness of p/s as well. □ 

4 An Application to Throughput Allocation 
on a Wireless Channel 

A central throughput allocation problem arising in the context of a wireless channel, such as 
in 3G technologies, is the following. There are n users 1, 2, ... n and the wireless router can be 
in any of m different states 1, 2, ... m whose probabilities, vr(j), can be estimated by sampling. 
Each user i derives utility at rate if it is connected to the router while the router is in state 
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j; the Ujj's are known. No matter what state the router is in, only one user can be connected to 
it. If user i is given connection for x^ < ir(j) of the time the router is in state j, for 1 < j < m, 
then the total utility derived by i is Vi = Sj=i u ij x ij- Clearly, we must ensure the constraint 
J27=i x ij — ^O); f° r each j. The question is to find a "fair" way of dividing the 7r(j)'s among 
the users. 

The method of choice in the networking community is to use Kelly's proportional fair scheme 
[Kel97], which entails maximizing J2i^°& v i subject to the constraints given above, i.e., solving 
the Eisenberg-Gale convex program [EG59]. Observe that the above setting can be viewed 
as a linear Fisher market with n users and m divisible goods. An elegant gradient descent 
algorithm for solving this convex program, given by David Tse [Tse] (see also [JPPOO]), was 
implemented by Qualcomm in their chip sets and is used by numerous 3G wireless basestations 
[And09]. However, this solution may at times allocate unacceptably low utility to certain users. 
This was countered by giving users the ability to put a lower bound on channel rates, say q 
for user i. This enhanced problem was solved by changing the objective function of the convex 
program to maximizing J2i^°s( v i ~ c «)> observe that this is precisely an instance of ADNB! 
However, now the gradient descent implementation ran into problems of instability, since it 
involved computing Uijj{v' i — q), where v • is the current estimate of vf, at intermediate points, 
the denominator may be too small or even negative. 

A different solution, proposed and implemented by researchers at Lucent [AQS05], was to 
introduce the constraints Vi > Ci in the Eisenberg-Gale program itself. The fairness guarantee 
achieved by this solution is unclear. Additionally, determining feasibility of the convex program 
now became a major issue [And09]. Instead, we have proposed experimenting with a heuristic 
adaptation of our combinatorial algorithm for ADNB, which will not have stability issues. 
As reported in the FOCS 2002 version of [DPSV08], an analogous heuristic adaptation of the 
DPSV algorithm was found to perform well on fairly large sized linear Fisher instances. 

5 Fisher's Model and its Extension via Flexible Budgets 

First we specify Fisher's market model for the case of linear utilities [BS00]. Consider a market 
consisting of a set of n buyers B = {1, 2, ... , n}, and a set of g divisible goods, G = {1, 2, . . . , g}; 
we may assume w.l.o.g. that there is a unit amount of each good. Let mi be the money 
possessed by buyer i, i £ B. Let Uij be the utility derived by buyer i on receiving one unit of 
good j. Thus, if is the amount of good j that buyer i gets, for 1 < j < g, then the total 
utility derived by i is 

9 

Vi(x) = ^ ' UijXij. 

i=i 

The problem is to find prices p = {pi,P2, ■ ■ ■ ,Pg} for the goods so that when each buyer is given 
her utility maximizing bundle of goods, the market clears, i.e., each good having a positive 
price is exactly sold, without there being any deficiency or surplus. Such prices are called 
market clearing prices or equilibrium prices. 
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The following is the Eisenberg-Gale convex program. Using the KKT conditions, one can 
show that its optimal solution is an equilibrium allocation for Fisher's linear market and the 
Lagrange variables corresponding to the inequalities give equilibrium prices for the goods (e.g., 
see Theorem 5.1 in [Vaz07]). 

maximize rrii log V{ (3) 

subject to Vi £ B : Vi = UijXij 

jeG 

Vz € B, Vj € G : Xij > 



Next, we introduce a flexible budget market as a modification of Fisher's linear case; this 
market will be used for solving ADNB. The utility functions of buyers are as before. The two 
main differences are that each buyer i now has a parameter a giving a strict lower bound on 
the amount of utility she wants to derive and buyers do not come to the market with a fixed 
amount of money, but instead the money they spend is a function of prices of goods in the 
following manner. Given prices p for the goods, define the maximum bang-per-buck of buyer i 
to be 



7,- = maxj < — - > . 
Now, buyer i's money is defined to be rrii = 1 + 

Let us say that set Si = argmax^ |^f| constitutes i's maximum bang-per-buck goods. Clearly, 
at prices p, any utility maximizing bundle of goods for i will consist of goods from Si costing 
rrii money. Again the problem is to find market clearing or equilibrium prices. Observe that 
in an equilibrium, if it exists, each buyer i will derive utility exceeding Cj. 



5.1 The Reduction 

An instance / of ADNB is transformed to a flexible budget market M. as stated above. 

Theorem 3 Instance I is feasible iff Ai is feasible. Moreover, if I and M. are both feasible, 
then allocations x and dual p are optimal for I iff' they are equilibrium allocations and prices 
for the flexible budget market M.. 



Proof : (=>) First assume that / is feasible and that allocations x and dual p are optimal 
for RNB game /. Then / must satisfy the KKT conditions for convex program (2). 
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By the second KKT condition, each good having a positive price is fully sold. Assume that 
yij > 0. Then, by the definition of 7, and the fourth KKT condition, 



u 



7i = — =Vi- Ci. 

Pi 

The money of buyer i at prices p in market M is defined to be m; = 1 + Cj/7^ The money 
spent by i in market A4 is: 

Xij Hij 

x ijPj = 2^ ~~ ; " 

jeG jeG '» 

1 ^ > 

= 2.^ X ij U ij = Z = 1 ^ Z = 1 ^ = m «- 

Furthermore, by the third and fourth KKT conditions, i buys only her maximum bang-per-buck 
objects, thereby getting an optimal bundle. This proves that x and p constitute equilibrium 
allocations and prices for market M. 

(-4=) Next, assume that A4 is feasible and that x and p are equilibrium allocations and prices 
for market Ai. Now, x is clearly feasible for program (2); we will show that x and p satisfy 
all the KKT conditions for this program. The first two conditions are obvious. 

Since i gets an optimal bundle of objects at prices p, 

Xij > => — = 7j. 
Pj 

Since i spends all her money, 

C{ ttik ^i 



mi = 1 + - = ^ x ijPj = ^ y ik — = -. 
Therefore, 7^ = — q. This gives the last two conditions as well. □ 



6 A Test for Equilibrium Prices 

We will present an efficient algorithm for solving an instance I of the game ADNB by first 
reducing it to a flexible budget market, say A4. In this section, we first give an efficient 
algorithm for the following simpler question: Given prices p = {pi, . . .p g } for the goods in A4, 
determine if these are equilibrium prices, and if so, find an equilibrium allocation. 

First, construct a directed network N(p) as follows. N(p) has a source s, a sink t, and vertex 
subsets B and G corresponding to the buyers and goods, respectively. For each good j £ G, 
there is an edge (s,j) of capacity pj, and for each buyer i € B, there is an edge (i,t) of 
capacity 777.3, where rrij = 1 + Cj/73 is i's money in M.. Recall that Si contains i's maximum 
bang-per-buck goods. The edges between G and B are precisely the maximum bang-per-buck 
edges, i.e., those (j, i) such that j G Si. Each of these edges has infinite capacity. 
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Lemma 4 Prices p are equilibrium prices for M. iff the two cuts (s, BuGUt) and (sUBL)G,t) 
are rain-cuts in network N(p). Moreover, ifp are equilibrium prices, then the set of equilibrium 
allocations corresponds exactly to max-flows in N(p). 

The proof of this lemma is straightforward using the transformation between a max-flow / in 
N(p) and an allocation x in M. given by Xij = f(j, i)/pj. The condition that (s, B U G U t) 
and (s U B U G, t) are min-cuts in network N(p), and hence saturated by /, corresponds to all 
goods being sold and all buyers' money being spent. The fact that (J, i) is an edge in N(p) 
iff j € Si ensures that buyers get only their maximum bang-per-buck goods. Clearly, one 
max-flow computation suffices to determine if prices p are equilibrium prices for M.. 

The next lemma gives the combinatorial object that yields equilibrium prices. Assume that 
p* are equilibrium prices, i.e., N(p*) satisfies the condition in Lemma 4. Let H be the unca- 
pacitated directed subgraph of N(p*) induced on B U G. 

Lemma 5 Given H, we can find p* in strongly polynomial time. 

Proof : Consider the connected components of H after ignoring directions on its edges. In 
each component, pick a good and assign it price p, say. The prices of the rest of the goods 
in this component can be obtained in terms of p. The bang-per-buck, and hence the money, 
of each buyer in this component can also be obtained in terms of p. Finally, by equating the 
money of all buyers in this component with the total value of all goods in this component, we 
can compute p. □ 

Given two ti-dimensional vectors with non-negative coordinates, p and q, we will say that p 
weakly dominates q if for each coordinate j, q,j < pj. The following characterization will be 
useful in Algorithm 3. 

Lemma 6 Assume that market A4 is feasible and p is its unique equilibrium price vector. Let 
q be a vector of positive prices such that (s, B U G U t) is a min-cut in network N(q). Then p 
weakly dominates q. 

Proof : Let us assume, for establishing a contradiction, that there are goods j such that 
qj > pj and yet (s,6UGUt) is a min-cut in N(q). Let 



Clearly, 9 > 1. 

Let T p and T q be the set of buyers who are interested in goods in S at prices p and q, respectively. 
Since S represents the set of goods whose prices increase by the largest factor in going from 
prices p to q, T q C T p . We claim that a buyer i G T q is not interested in any goods in G — S at 
prices p (because otherwise at prices q, i will not be interested in any goods in S, since their 



6 = max < — 
jeG [pj 




and S = {j € G \ qj = Opj}. 
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prices increased the most). Therefore, in any max-flow in N(p), all flow going through nodes 
in T q must also have used nodes in S. Therefore, 



jeS ieT q v '■ 

where j t is the maximum bang-per-buck of buyer i w.r.t. prices p. Multiplying this inequality 
by 6 and using the fact that 9 > 1 we get, 

jeS ieT q v '* 7 i& q v '» 

Observe that the maximum bang-per-buck of buyer i G T q w.r.t. prices g is 7j/0. Therefore, 
the last inequality implies that w.r.t. prices g, the total value of goods in S is strictly more 
than the total value of money possessed by buyers in T q . On the other hand, since in N{q) all 
flow using nodes of T q goes through S, we get that (s, B U G U t) is not a min-cut in network 
N(q), leading to a contradiction. □ 

Next assume that Ai is an arbitrary flexible budget market, not necessarily feasible. We will 
say that prices q are small if q is a positive vector and (s, B U G U t) is a min-cut in network 
N(q). Observe that in this case, each good j must have an edge (J, i), for some buyer i, incident 
at it. By Lemma 6, if M is feasible and prices q are small, then they are weakly dominated 
by the equilibrium prices, p. Observe however that the contrapositive of Lemma 6 does not 
hold, i.e., p may dominate positive prices q, yet (s, BuGUt) may not be a min-cut in network 
N(q). 

The background given so far will suffice to read Section A, which gives an algorithm that 
converges to the solution of a given feasible instance of ADNB in the limit. This section may 
serve as a suitable warm-up, since our polynomial time algorithm is quite involved. 



7 Determining Feasibility 

In this section, we will address the question of determining whether the given flexible budget 
market, M, is feasible. We will give a characterization of feasible markets and we will derive 
conditions that yield a proof of infeasibility. 

Clearly, if we can find small prices p and a max-flow / in network N(p) such that the flow 
gives each buyer i strictly more than q utility, then M. is feasible. We first show, using the 
notion of balanced flows, that this test of feasibility is in fact a property of prices p only. 



7.1 Balanced flows 

We will follow the exposition in [Vaz07] and refer the reader to this chapter for all facts stated 
below without proof. For simplicity, denote the current network, N(p), by simply N. Given a 
feasible flow / in N, let R(f) denote the residual graph w.r.t. /. Define the surplus of buyer 
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i w.r.t. flow / in network N, 9i(N, /), to be the residual capacity of the edge (i, t) w.r.t. flow 
/ in network N, i.e., rrii minus the flow sent through the edge (i,t). The surplus vector w.r.t. 
flow /is defined to be 9(N,f) := (Oi(N, f),9 2 (N, /),..., n (N, /)). Let ||v|| denote the l 2 norm 
of vector v. A balanced flow in network iV is a flow that minimizes ||6>(A, /)||. A balanced flow 
must be a max-flow in N because augmenting a given flow can only lead to a decrease in the 
I2 norm of the surplus vector. 

A balanced flow in N can be computed using at most n max-flow computations. It is easy to 
see that all balanced flows in N have the same surplus vector. Hence, for each buyer i, we can 
define 9i(N) to be the surplus of i w.r.t. any balanced flow in N; we will shorten this to 9i 
when the network is understood. The key property of a balanced flow that our algorithm will 
rely on is that a maximum flow / in N is balanced iff it satisfies Property 1: 

Property 1: For any two buyers i and j, if 9i(N,f) < 9j(N,f) then there is no path from 
node i to node j in R(f) — {s, t}. 

Balanced flows play a crucial role in both stages of our algorithm; moreover, they have multiple 
uses. In Section 8.4, after stating the full algorithm, we state the various uses of this notion. 

7.2 A characterization of feasibility 

Let p be small prices and let (9±, . . . , 9 n ) be the surplus vector of a balanced flow in N(p). We 
will say that p are feasible prices if for each buyer i, 9i < 1. 

Lemma 7 Market M. is feasible iff there are feasible prices for it. 

Proof : If M. is feasible, its equilibrium prices are feasible, since for each buyer i, 9i = 0. 
Next, assume that p are feasible prices for M. By definition, the flow sent on edge (i,t) in a 
balanced flow in N(p) is 

mi — 9i > rrii - 1 = ( 1 + — j — 1 = — . 

V 7* / li 

The utility accrued by i from this allocation is 7j(mj — 9i) > c«. Hence M is feasible. □ 

Observe that if a max-flow / in network N(p) gives each buyer i strictly more than Cj utility, 
then so will a balanced flow in N(p). Hence, feasibility is a property of prices p only. 

Rather than working with 9{, it will sometimes be more convenient to work with 0j — 1. Hence, 
w.r.t. small prices p, let us define the 1-surplus of buyer i to be /3j = 9i — 1, where (#1, . . . , 9 n ) 
is the surplus vector of a balanced flow in N(p). Now, another definition of feasible prices is 
that they be small and for each buyer i, ^ i < 0. 

7.3 Proof of infeasibility 

We will provide two ways of establishing infeasibility of the given market. The first is via the 
dual of an LP whose optimal solution tells us if M. is feasible and the second is via the dual of 
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convex program (2). Given an infeasible market, Stage I of our algorithm will terminate in a 
way that yields both proofs. 

The game is feasible iff there is a point v G J\f such that for each agent i G B, Vi > Ci. In order 
to capture feasibility via a linear program, let us restate as follows: the game is feasible iff 

max min : (v j — Cj) > 0. 
veAf i£ B 

Observe that the expression on the left hand side is the optimal objective function value of 
LP (4). Hence, the game is feasible iff the optimal solution to LP (4) is greater than zero. 
Clearly, this LP is maximizing t; however, in order to obtain a convenient dual, we will write 
it as minimizing —t: 

minimize —t (4) 
subject to Vi G B : ^ UijXij > Cj + t 

Vj G G : -]T^>-1 

ieB 

Vi G B, Vj G G : Xij > 

Let yiS and Zj's be the dual variables corresponding to the first and second set of inequalities, 
respectively. The dual program is: 

maximize ciyi — V] Zj (5) 

iei jeG 

subject to \/i £ B, Vj G G : n^yj — Zj < 

VieB : yi >0 
Vj G G : z,- > 

Lemma 8 // there exist prices p s.t. J2jecPj > a^d X^eB ^ 0, i/ien the given game is 
infeasible. 



Proof : W.r.t. prices p, compute the maximum bang-per-buck, j i , of each buyer i. By 
definition of maximum bang-per-buck, 

Vi G B, Vj G G : 7i > ^. 
Let n = J2ieB l hi- Since T,jeGPj > °> M > °- 

Next, consider prices q, where for each i G G, (fy- = Since all the prices have been scaled 

by the same factor, the network remains unchanged. Clearly, the maximum bang-per-buck of 
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buyer i w.r.t. q is 7^ = and 

V* E B, Vj € G : ^ > ^- 

Let i/i = 1/7-, for i G B, and £j = (/j, for j G G. We will show that (y, z) is a feasible solution 
for the dual LP (5). The first set of inequalities is established by noting that 

Mi G B, Vj E G : 7- > — hence UijUi < Zj. 
Next we show that the equality constraint holds: 

Let a- = Cj/7^ and let j3[ be the 1-surplus of buyer i w.r.t. prices q. The objective function 
value of the dual solution (y, z) is 



E ^ ~ E z 3 = E a 'i ~ E «i = E # = 

ieB jeG i&B jeG i£B 




Therefore, at optimality, —t > 0, i.e., t < 0, hence establishing infeasibility of the game. □ 

Next, we derive a condition under which the following convex program, which is the dual of 
(2), has an unbounded solution, hence proving infeasibility of the primal. This dual program 
was given by Devanur [DevlO]. Its variables are qj's and j/j's. 

minimize E ^ ~ E ^ ~ E lo g(^) ( 6 ) 

jeG i€B i<=B 

subject to Mi G B, Vj E G : #j > u^-yj 

Lemma 9 If there exist positive prices p such that network N(p) can be partitioned into 
two, induced on (B',G f ) and ((B - B'),(G - G')), where B' C B and G' C G, such that 
J2ie(B-B') /3j > and Mi G (B — B'), Vj G G' : Uij = 0, then the given game is infeasible. 

Proof : Observe that 

E ft = ( E *i) - |B - J^l = ( E rn)-( £ p^-IS-B'l 

iG(B-B') \ie(B-B>) J \ie(B~B>) J \je(G-G') J 

= ( £ ? )-f E ») 

\*e(B-B') 'V Vje(G-G') / 

It is easy to see that setting q.,' to pj and to 1/7$ gives a feasible solution to program (6). 
Multiply the prices of all goods in (G — G') by x and let x — > 00. Observe that because of the 
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condition Vi G (B — B r ), Vj G G' : = 0, no new edges will be introduced in the network. 
Hence, the updated setting of q^s and y^'s still yields a feasible solution to the dual. 

As x — >■ 0, for each i £ (B — B'), log(7j — > — oo. Also, ^je(G-G') 9j ~~ Sie(B-.B') c «y« is either 
or tends to — oo. Hence the entire objective function tends to — oo. Therefore, the dual is 
unbounded and hence the primal is infeasible. □ 



8 Details of the Algorithm for ADNB 

We will impose the following condition throughout; by Lemma 6, it will ensure that prices are 
always small 2 . 

Invariant: W.r.t. current prices, p, (s, B U G U t) is a min-cut in network N(p). 

It is easy to see that the prices found by Initialization satisfy the Invariant. 

Let / be a balanced flow in N(p). Since the Invariant is always maintained, for each buyer i, 
9i > and hence /3j > — 1. In the algorithm, we will change prices of a well-chosen set J of 
goods as follows. Multiply the price of each good in J by a variable x and initialize z to 1. In 
Stage I, we will decrease x and in Stage II we will raise x until the next event happens. 

In the next lemma, we will assume that J = G and we will study how the 1-surplus of buyers 
changes as a function of x. Define x ■ f to be the flow obtained by multiplying by x the 
flow on each edge w.r.t. /. Let denote i's 1-surplus w.r.t. flow x ■ p. Let B 1 be the 

set of buyers having negative 1-surplus w.r.t. prices p. If B' = 0, define b = oo else define 
b = minj e s'{— l//3j}. Observe that in both cases, b > 1. 

Lemma 10 Flow x ■ f is a balanced flow in N(xp) for < x < b, and for each i G B, Pi(x) = 



Proof : Since the Invariant holds and / is a max-flow in N(p), the cut (s, J U I U t) is 
saturated by /, and hence by x ■ f in N(xp). Next we show that x ■ f is a feasible flow in 
N(xp), i.e., for each buyer i G B, edge (i,t) is not over saturated. Now, /3j = ccj — f(i,t). 
Therefore, the surplus on edge (i,t) w.r.t. flow x ■ f is 1 + x(cti — f{i,t)) = 1 + x^ i > for 
< x < b. Hence, edge (i,t) is not over saturated. Furthermore, /3j(x) = xf3 i . Finally, since / 
satisfies Property 1 in N(p), x ■ f satisfies it in N(xp), thereby showing that it is a balanced 
flow. □ 



8.1 Details of Stage I 



Algorithm 1 gives the pseudo code for Stage I. In this section, we give the subroutines used by 
this stage and Section 8.3 gives formal definitions of the predicates used in the While loops. A 

2 The following power point presentations may make the algorithm easier to understand: 
http:/ /www. cc.gatech.edu/ vazirani/Waterlool.ppt 
http://www.cc.gatech.edu/ vazirani/Waterloo2.ppt 
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run of Stage I is partitioned into phases, which are further partitioned into iterations. In Stage 
I, an iteration ends when a new edge is added to the network. A phase ends either when the 
condition of Step 7 holds or if for some i G I, fa > 0. In each iteration, the algorithm computes 
a balanced flow in the current network, N(p). 

We establish the following notation. For J C G, define p(J) = J2jt=jPj an d T(J) = {i G 
B | 3j G J si. (j, i) G N(p)}. Similarly, for I C B, define a(7) = X^eJ a i> m(J) = J2iei m i 
and T(I) = {j £G \3i e B s.t. (j,i) £ N(p)}. 

The sets -B c and G c denote the current sets of buyers and goods being considered by the 
algorithm. These sets are initialized to B and G, respectively. At any point the algorithm, 
B = B c U B' and G = G c U G', where B' and G' are the sets of adaptable buyers and goods, 
respectively; their purpose is explained below. B' and G' are both initialized to 0. As the 
algorithm proceeds, buyers are moved from B c to B' and goods are moved from G c to G' . 

The subroutines used in Stage I are: 

• Find sets (I): Sets I C B c and J C G c are initialized as follows. 

I<- argmin{/?J and J <- (r(J) - r(B c - /)). 

ieB c 

Observe that J consists of goods that are the maximum bang-per-buck goods of buyers 
in I only. All edges from goods in G c — J to buyers in I are removed; this is justified in 
Lemma 15 below. 

• Update sets(I): Find the set, I', of all buyers in B c — I such that there is a residual 
path from a buyer in / to a buyer in I'. Update 

I-G-(IUI') and J 4- {T(I) -T(B C - I)). 

All edges from goods in G c — J to buyers in I are removed. Once again, this is justified 
in Lemma 15. 

Assume that (j,i), j G J, % G (B c — I) is the new edge added to Si in the current iteration. 
Observe that if I' = 0, then all the flow from j, which was going to buyers in I before the 
addition of this edge, must go to i, since there is no residual path from / to i. Accordingly, 
Update sets (I) will move good j from J to G c — J. As soon as the prices of goods in J are 
reduced by an infinitesimally small amount (by decreasing x) buyers in I will not be interested 
in good j anymore. 

Lemma 11 In Stage I, at the start of each iteration, for each buyer i G I there is a good j G J 
such that edge ( j, i) is in the network. 



Proof : Since fa < 0, the balanced flow must be sending flow on some edge (j, i). If an edge 
(j, i') to a buyer i' G (B c — I) is also present in the network, then there will be a residual path 
from i to i! , violating Property 1. Therefore, there is no such edge and j G (r(J) — T(B C — I)), 
proving the lemma. □ 
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We now explain the purpose of the sets B' and G' . Once a good is moved into G', its price 
gets frozen until the end of Stage I. At any point in Stage I, these sets satisfy the following 
properties: 

1. W.r.t. the frozen prices of goods in G', for each buyer i G B', (3 i < 0. 

2. Buyers in B are totally uninterested in goods in G' at any price, i.e., for every % G B and 
every j G G', Uij = 0. Hence, as prices of goods in G' are decreased, no edge from B to 
G' will ever enter the network. 

The reason for the name "adaptable" is that as far as determining feasibility or infeasibility 
goes, buyers in B' can be made consistent with the outcome of the remaining buyers. Thus, 
if Vi G B c , /3j < 0, by assigning the frozen goods in G' their prices at the time of freezing, we 
can ensure that Vz G B, f3 i < 0. The resulting price vector is clearly feasible. In Step 9 in 
Algorithm 1 we have refer to this process as "restoring prices of adjustable goods." 

If on the other hand J2ieB c @i — 0> then we can give a proof of infeasibility in one of two ways. 
First, by lowering the prices of all goods in G' to zero, which can be done without introducing 
any new edges in the network, we can ensure that Vi G B', fi i = 0, thereby ensuring that these 
buyers don't affect the sum of /3j's. The conditions of Lemma 8 now hold and yield a proof of 
infeasibility. Second, by assigning the frozen goods in G' their prices at the time of freezing, 
we can ensure that the prices of all goods are positive and the conditions of Lemma 9 now 
hold to yield a different proof of infeasibility. 

At the start of a phase, the set / C B c of buyers having the smallest f3 values is identified. 
The goods they desire are put in set J. If at any point, I and J are found to be adaptable, 
the algorithm updates B' and G 1 and the phase comes to an end. Otherwise, the algorithm 
lowers the prices of goods in J until a new edge (j, i), with j £ J and i G (B c — I) is added to 
the network. On recomputing a balanced flow, either f3 i becomes negative, if so i moves into I 
and the iteration comes to an end, or for some buyer(s) %' G I, (3 V increases. If for some buyer 
i G /, /3j becomes non-negative, the phase comes to an end. Lemma 12 shows that eventually, 
either all buyers are rendered good or the conditions of Lemma 8 start holding. In the former 
case, the frozen prices of goods in G' are restored and the algorithm moves on to Stage II to 
find equilibrium prices. 

One way to view the operation of Stage I is as a tug-of-war between two sets of buyers: the 
good buyers and the rest. The algorithm decreases the prices of goods desired by buyers in 
/, thereby increasing their /3j's. This helps towards reaching the infeasibility condition stated 
above. However, as new edges enter the network and a balanced flow is recomputed, buyers 
may move between the 2 sets. 

Observe that in each iteration, the algorithm needs to compute the largest value of x at which 
a new edge is added to the network. For any one edge this is straightforward; taking the 
maximum over all relevant edges gives the required value. 

Lemma 12 Stage I must terminate with either a feasible price vector or a proof of infeasibility. 

Proof : By Lemma 25, Stage I must terminate. If at this point, Vi G B : (3^ < 0, a feasible 
price vector has been found. Otherwise, J2i<=B c Pi — and hence 3i G B c : p i > 0. This 
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buyer must be in B c — I, since all buyers in I satisfy < 0, and the goods this buyer desires 
must have positive prices. Now, the conditions of Lemma 8 can be made to hold by setting 
the prices of goods in G' to zero. □ 



8.2 Details of Stage II 

Algorithm 2 gives the pseudo code for Stage II. In this section, we give the subroutines used 
by this stage and Section 8.3 gives formal definitions of the predicates used in the While 
loops. If Stage I terminates with a feasible price vector p, the algorithm moves to Stage II 
to find equilibrium prices. Since p is small, Stage II needs to raise prices of goods to get to 
the equilibrium. Another way to view the situation is that since buyers have surplus money, 
Stage II needs to change prices in such a way that the surplus drops to zero. Are both 
these requirements compatible, i.e., will the surplus of buyers decrease by raising prices? The 
following lemma clarifies this crucial point in the simplified setting of Lemma 10, i.e., prices of 
all goods are raised; of course, Stage II will raise the prices of well-chosen subsets of G. 

Lemma 13 If in the setting of Lemma 10, prices p are feasible and x is raised without violating 
the Invariant, then the surplus of each buyer decreases and the resulting price vector is still 
feasible. 

Proof : Since p is feasible, for each buyer i, f3 i < 0. Clearly, if x > 1, x ■ f3 i < j3, h i.e. the 
surplus of buyer i decreases. Moreover, the property that the f3 of each buyer is negative is 
preserved. Hence the resulting price vector is still feasible. □ 

A run of Stage II is also partitioned into phases, which are further partitioned into iterations. 
An iteration ends when a new edge is added to the network and a phase ends when a new 
set goes tight. We will say that S C G is a tight set if the total price of goods in S exactly 
equals the money possessed by buyers who are interested in goods in S, i.e., p(S) = m(T(S)). 
Clearly, if S is tight, buyers in T(S) must have zero surplus and hence have /3j = —1. In each 
iteration, the algorithm computes a balanced flow in the current network, N(p). 

The subroutines used in Stage II are: 

• Find sets(II): Sets I C B and J C G are initialized as follows. 

I <- argmax{#j} and J T(I). 

All edges are removed from goods in J to buyers in B — 7; this is justified in Lemma 15 
below. 

• Update sets(II): Find the set, I', of all buyers in B — I that have residual paths to 
buyers in /. Update 

I <- (III I') and J«-r(I). 
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All edges are removed from goods in J to buyers in B — I. Once again, this is justified 
in Lemma 15. 



Observe that if is the new edge added to Si, then good j must move from G — J to J, 
whether or not I' = 0. The choice of set J above ensures that if the prices of goods in J are 
increased by an infinitesimally small amount (by increasing x as stated in Algorithm 2), there 
is no change in the maximum bang-per-buck goods of buyers in B. 

Lemma 14 In Stage II, at the start of each iteration, for each buyer i G (B — I) there is a 
good j G (G — J) such that edge (j, i) is in the network. 

Proof : Since B\ < 1, the balanced flow must be sending flow on some edge (j, i). If j G J, 
then there will be a residual path from i to a buyer in /, violating Property 1. Therefore, 



Lemma 15 In Stage I (Stage II), the Invariant holds after all edges from goods in G — J (J) 
to buyers in I (B — I) are removed. 

Proof : The idea of the proof is the same for both statements. In Stage I, right after 
Update sets (I) is executed, there are no residual paths from I to B — I. Therefore, by 
Property 1, any edges from G — J to I could not be carrying any flow and hence their removal 
will not affect the Invariant. 

In Stage II, right after Update sets (II) is executed, there are no residual paths from B — I 
to /. Therefore, by Property 1, any edges from J to B — I could not be carrying any flow and 
hence their removal will not affect the Invariant. □ 

In each iteration, we need to compute the smallest value of x at which a new edge is added to 
the network or a new set goes tight. The former computation is the same as in Stage I. Let 

the smallest value of x at which a new set goes tight be x* . Let b = minj e / j — j . Clearly, 

b > 1. Using Lemma 10, proved in Section 7, for x in the range 1 < x < b, we prove below 
that x* = b. 

Lemma 16 x* = b. 

Proof : By definition of b, for 1 < x < b, for each i €. I, the surplus of i will be 1 + xfli > 0, 
since xf3 i > —1. Thus, each edge (i,t) will have positive surplus, implying that there are no 
tight sets. 

Next, assume that x = b. Let 



je(G-J). 



□ 



T = 




— = b \ and S = { j G J | /(j, i) > 0, for some i G T}. 



i 
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Since / is a balanced flow in N(p), there cannot be an edge (j, i) for j G S and i G (7 — T) in 
N(p). This is so because otherwise there would be a path from T to i in the residual graph, 
contradicting Property 1 (observe that fi i < —1/6). Therefore, T(S) = T. Moreover, for i G T, 
the surplus on edge (i,t) w.r.t. flow x • / in N(xp) is 1 + x(— 1/b) = 0. Hence S is a tight set 
in network N{xp) for x = b. □ 



8.3 Predicates used in While loops 

In Step 2 (Stage I), "a proof of feasibility is reached" when Mi G B, /3 f < 0, and "a proof of 
infeasibility is reached" when X^eB A — 0- Thus "a proof of feasibility or infeasibility is not 
reached" is satisfied iff -i((£ ieB ft > 0) V (Vt eB, ft< 0)). 

In Step 4 (Stage I), "5 - 7 desire J" is satisfied iff 3i € (73 - 7), 3j e J : u^- > 0, and 
"buyers in 7 have small surplus" is satisfied iff Vz G 7, /3j < 0. 

In Step 1, "a buyer in B has surplus money" is satisfied iff 3i G 7?, 0j > 0). 

In Step 3, "no set in J is tight" is satisfied iff ~^(3S, C S C J si. 5 is tight). 

8.4 The role of balanced flow 

Besides being used for defining the central notion of feasible prices, balanced flow plays the 
following three, rather diverse, crucial roles in both stages of our algorithm. 

1. Ensure that edges, that need to be removed as prices of goods in J are raised, did not 
carry any flow and hence their removal would not violate the Invariant; this is argued in 
Lemma 15 

2. Ensure that in each iteration, buyers entering 7 in Stage I (Stage II) have sufficiently 
large this is established in Lemma 21 (Lemma 31). 

3. Prove that sufficient progress is made in an iteration and hence in a phase. This is 
established in Lemma 22 for Stage I and Lemma 32 for Stage II. 

As stated in the Introduction, balanced flow could have been defined without resorting to the 
I2 norm - as a max-flow that makes the surplus vector lexicographically smallest, after its 
components are sorted in decreasing order (and hence making the components as balanced 
as possible). It is easy to prove Property 1 with this definition as well. The first two roles 
listed above make use of Property 1 only. On the other hand, the third role uses the definition 
of balanced flow via the 1% norm and as argued in Section 11, the use of the I2 norm seems 
indispensable. 
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Algorithm 1 (Initialization and Stage I of the Algorithm 


for ADNB) 


1. 


Initialization: 






(i) Vi G -B : mi <- 1. 






(ii) Use the DPSV algorithm to compute equilibrium prices, p 






(hi) Vi£B: ^ 1 + ^. 

I i 






(iv) B c ^B; G c ^ G. 






(v) B' ^0; G'<-0. 






(vi) Compute a balanced flow in N(p). 






Stage I 




2. 


fiVeio Phase) While a proof of feasibility or infeasibility is not reached do: 


i. 


r inQ sets^i^. 




4. 


(New Iteration) While B c — I desire J and buyers in / have 


small surplus do: 


0. 


Multiply the prices of goods in J and as of buyers in / by x. 






Initialize x <— 1, and decrease x continuously until: 






A new edge (j, i) enters Si, for j G J and i G (B c — I). 






Add (j, i) to iV(p) and compute a balanced flow in it. 






Update sets(I). 




p. 

D. 


End (End Iteration) 




"7 


If Mi G (B c - I), Vj G J : u^- = 0, then: 






Declare / and J adaptable, i.e., 






G'^(G'UJ) and G c <- (G c - J). 






B'^(B'UI) and B C ^{B C -I). 




8. 


End (£nd P/iase; 




9. 


If Vi G Be, ft < 0, then: 






Restore prices of adaptable goods in G' . 






Compute a balanced flow in N(p). 






Go to Step 1 in Stage II. 




10. 


Else (i.e., J2ieB c Pi ^ 0)- output "The game is infeasible". 






HALT. 
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Algorithm 2 (Stage II of the Algorithm for ADNB) 

1. (New Phase) While a buyer in B has surplus money do: 

2. Find sets(II). 

3. (New Iteration) While no set in J is tight do: 

4. Multiply prices of goods in J and as of buyers in I by x. 

Initialize x 4— 1, and raise x continuously until: 

A new edge enters Si, for j G (G — J) and % G I. 

If so, add (j, i) to N(p) and compute a balanced flow in it. 
Update sets(II). 

5. End (End Iteration) 

6. End (End Phase) 

7. Output the current allocations and prices. 
HALT. 



9 Running Time Analysis 

We first define some parameters of the given problem instance. Recall that g = \G\ and 
n = \B\. Let U = maxi£B,jeG{ u ij}, C = maxjgsQ, and A = nCU n . Observe that program 
(2) with all Cj = is the same as the convex program for a linear Fisher market with all buyers 
having unit money. Hence, Theorem 2 gives a lower bound on the price of a good computed 
in Initialization. Let this lower bound be l/ji, fi G Z + . 

The following enhanced version of Lemma 10 will be needed in both stages. 

Lemma 17 Let f be a balanced flow in network N(p). Then, for < x < b, the flow x ■ f is 
a balanced flow in N{xp). 

Proof : For i,j G B assume that 1 + xf3 { < 1 + x(3j. Since x > 0, 1 + t < 1 + (3j, i.e., 
w.r.t. flow / in N(p), the surplus of i is smaller than that of j. Since / is a balanced flow in 
N(p), by Property 1, there is no path from i to j in the residual graph. Therefore, w.r.t. flow 
x ■ f in N(xp) also there is no path from i to j in the residual graph. Therefore, flow x ■ f in 
N(xp) satisfies Property 1 and hence is a balanced flow. □ 

9.1 Stage I 

Throughout Stage I, we will consider a partitioning of B c into two sets, B\ and B2, containing 
buyers having /3 i < and f3 { > 0, respectively. For Stage I, we will work with the following 



23 



potential function: 



As Stage I proceeds, buyers move from £> c to 5', and within B c between the sets B\ and B2. 
For this reason, it will be convenient to define $ using an n-dimensional vector, ijj, called the 
associated vector of network N. The i-th component of this vector, ipi, is P i for i G B±, and is 

for i E (B' U B2). Hence, an alternative definition for the potential function is: 

Lemma 18 In Stage I, a phase consists of at most ng iterations. 

Proof : Observe that if in Update sets(I), I' = 0, then a good must move from J to G c — J. 
Otherwise, a buyer must move from B c — I to /. Clearly, there can be at most \J\ < \G C \ 
contiguous iterations of the first type and a total of at most \B C — Iq\ < \B C \ iterations of the 
second type, where Io is the set / at the start of the phase. □ 

The central fact established below is that $ drops by a factor of (1 — l/(gn 2 )) in a phase 
(Lemma 23). Towards this end, assume that a given phase consists of k iterations. Let Io 
denote set / at the start of the phase and let // denote the set I at the end of the l-th iteration, 

1 < I < k. Assume that at the start of this phase, maxjg^d/^l} = 6 = 5q. Let 

Si = min{|/3j|}, for 1 < I < k, and 5j~ = 0. 

ieii 

As we will see in this section, the potential function <3? drops monotonically in each iteration 
in the phase. Within an iteration, we will account for the drop in two steps. First, as prices of 
goods in J are reduced, by Lemma 10 the /3j's of buyers i E I increase, leading to a reduction 
in Second, when a new edge (j, i), with j G J and i € (B c — I), is added to the network, 
the flow becomes more balanced, leading to a further drop. We will account for these two 
reductions separately, via different arguments (see Lemma 22). For the first step, we work 
with the l\ norm, establishing an increase in X^eBi Pi- m the second step, J2ieB 1 Pi wm n °t 
change if i G (B\ — I). Instead, we establish a decrease in ||^|| 2 using an I2 norm based 
argument. We observe that the latter argument is difficult to apply to the first step since the 
money of buyers changes as prices change. Also, we do not know of a simple one step argument 
that accounts for the entire reduction in an iteration. 

Next, we prove a key fact that accounts for the second decrease. Just before new edge (j, i) 
is added to Si, let N be the network and p be the prices of goods. Let N' be the network 
obtained by adding this edge to N; of course, the prices remain unchanged. Let / and /* be 
balanced flows in iV and N', respectively, and let tpi and be their associated vectors. 

Lemma 19 ||V>|| 2 - ||V>1| 2 > EheB, tyh ~ ■ 

Proof : Since the Invariant holds and the prices are unchanged, / and /* have the same 
value. Therefore, flow f* — f will consist of circulations. Since / is a balanced flow, all these 
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circulations must use the edge (j, i), because otherwise a circulation not using edge could 
be used for making / more balanced. These circulations will have the effect of increasing 
the surplus of certain buyers in /, say i\, for 1 < I < k, and decreasing the surplus of buyer 
ie{B c -I). Let fa - fa = 5, and for 1 < I < k, /?* - fa { = S t . Then, £f=i 5 t = 5. 

For each buyer ii, 1 < i\ < k, there is a path from i\ to i in the corresponding circulation and 
hence there is a path from i to i\ in the residual graph w.r.t. flow /*. Since /* is balanced, by 
Property 1, the surplus of buyer i is at least as large as that of i\. Therefore, fa > /?* . 

In going from N to N', the iph values can change only for h = i, and h = ii, for 1 < I < k. We 
will consider 3 cases. 

Case 1: fa > 0, i.e., i € B 2 , and fa > 0. 

In this case, ipi = ip* = and the lemma is obvious. 

Case 2: fa > 0, i.e., i e B 2 , and fa < 0. 

Let a = — fa. In this case, ipi = and ip* = —a. 

Clearly, a < £f=i <S, Since fa >fa t ,a<bi-5i. Now, 

IHI 2 - lllHI 2 = fo 2 + E^-fa 2 + E(^-^) 2 )=-a 2 + E(2fc/-W 

V 1=1 ) \ 1=1 ) 1=1 

> -a 2 + E (2a + <S,)*l > -a 2 + E <*? + 2o E ^ ^ E 

z=i z=i i=i i=i 

where the first inequality follows from a < bi — 5i and the third one follows from a < Yli=i 
Case 3: fa < 0, i.e., i € {B\ - I). 

Clearly, in this case, fa < 0. Substitute a = —fa, and for 1 < i\ < k, substitute bi = —fa r 
Now, by Lemma 20, to get \\ip\\ 2 - \\ip*\\ 2 > 5 2 . Clearly, 5 2 > £*= 

l $i ) giving the lemma. □ 

Lemma 20 Let 5, 5 t > 0, I = 1, 2, . . . , k, with 5 = £?=i<*J- If a + S < b t - S h for I = 1,2, ... ,k 
then 

\\{aM,b 2 ,...,b k )\\ 2 -\\{a + 5,^-5^2 -5 2 ,...,b k -5 k )\\ 2 >5 2 . 
Proof : 

(« 2 + E^)-((« + ^ + E(^-^) 2 ) 

> ((a + 5 - 5) 2 + E(6i " ^ + <^) 2 ) " ((a + <5) 2 + E (&i " *i) 2 ) 
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□ 

Let tp° denote vector ib at the start of the phase and ip l denote tp at the end of iteration I, for 
1 < I < k. 

Lemma 21 In the l-th iteration, there is a buyer i € such that decreases by at least 
(<$,_i - <$,), for 1 < I < k. 

Proof : By the definition of set I 1 in procedure Update sets (I) and Property 1, there is 
a buyer i £ which achieves minj G / ; {|/3j|} at the end of iteration /. Clearly, ji i increases 
(and hence decreases) by at least (<5/_i — Si) in the Z-th iteration. □ 

Lemma 22 For 1 < I < k, 

ll^f-II^II^W-l-*!) 2 . 

Proof : We first prove the statement for 1 < I < k. By Lemma 21, there is a buyer i G 
such that P i increases by at least — Si) in the Z-th iteration. Let us split this increase into 
two parts, the increase due to decrease in the prices of goods in J and that due to a new edge 
entering the network. Let these be a and b, respectively. Therefore, a + b = <^_i — S^ 

Let ip' be the vector ip just before the new edge is added to the network in iteration I, i.e., 
right after all the decrease in prices of J has happened. As prices in J decrease, the beta's of 
buyers in / increase, each leading to a decrease in H^'H 2 ; clearly, the beta's of buyers in B c — I 
remain unchanged. Let c be the value of beta of buyer i at the beginning of iteration /. Then, 

\\^- 1 f-H'\\ 2 >c 2 -(c + a) 2 = a 2 -2ac. 

By Lemma 19, 

||^/||2_||^||2> 6 2_ 

Adding the two we get 

||2 _ ||^||2 > a 2 _ 2ac + b 2 > ( y + y/ )2 > _ ( j j)2j 

where the second last inequality follows from the observation that b < —c. 

Finally, in the k-th iteration, there is a buyer i £ Ik-i whose ipi changes from < to 0. 
Therefore, 

||^|| 2 > yis^-Sk) 2 , 

since S^-i < —fli and Sk = 0. □ 
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Lemma 23 In a phase in Stage I, the potential drops by a factor of 




Proof : Now, \\ip°\\ 2 — HV^II 2 can be written as a telescoping sum of k terms, each of which is 
the decrease in the potential in one of the k iterations. Lemma 22 gives a lower bound on each 
of these terms. The total lower bound is minimized when each of the differences (<5/_i — Si) is 
equal. Now using the fact that So = S and Sk = 0, we get: 

u \?-u k \\ 2 > 5 ^- 

Finally, since \\i^> \\ 2 < nS 2 , and by Lemma 18 k < ng, we get: 

□ 

Lemma 24 At any point in Stage I, if J2i^B 2 Pi > ®> then 

Proof : First observe that the maximum bang-per-buck of buyers in B2 remains unchanged 
throughout Stage I, and is determined by prices of goods found in the Initialization. Let 
J2 = r(£?2)- By Property 1, there is no flow from a good in J2 to a buyer in Bi, and therefore, 
all flow from J 2 must go to buyers in B 2 . Therefore, 

9 t -\B 2 \ = j2^-T,pj- 

i€B 2 ieB 2 ieB 2 j€J 2 

Now if 2~ZieB 2 Pi > 0) then the denominator of this sum is a product of at most n u^s and at 
most g p/s, and is therefore bounded by U n /j, 9 , proving the lemma. □ 

Lemma 25 The execution of Stage I requires at most 

O (n A g 2 {n log U + g log 
max-flow computations. 

Proof : By Lemma 23, the square of the potential drops by a factor of two after 0(n 2 g) 
phases. At the start of the algorithm, the potential is at most n. 
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If at any point, J2ieB 2 Pi = 0> i- e -' G -B2 : Pi = 0, each subsequent phase must end up in 
Step 5 and not Step 6, i.e., some buyers will be removed from consideration. Therefore, from 
this point on, at most n more phases are needed for Stage I to terminate. 

Next, let's assume that X^eB 2 Pi throughout Stage I. If so, by Lemma 24, once the potential 
drops below 1 - 2 , the phase must end (in Step 7). Therefore the number of phases is 

0(nVog(£/ 2 V 9 ))- 

By Lemma 18 each phase consists of at most ng iterations and each iteration requires n max- 
flow computations for finding a balanced flow. The lemma follows. □ 

9.2 Stage II 

When Vi G B : P i < 0, the algorithm starts with Stage II. Since in this stage the algorithm 
only raises prices of goods (i.e., increases x), by Lemma 13, Vi G B : P i < holds until 
termination. 

In this section, we will work with the O^s of buyers, rather than their /3j's. Thus, throughout 
Stage II, Vi G B : 0, < 1. For Stage II, we will work with the following potential function: 

ieB 

Lemma 26 In Stage II, at the termination of a phase, the prices of goods in the newly tight 
set must be rational numbers with denominator < A. 

Proof : Let S be the newly tight set. Consider the subgraph of the network induced on the 
bipartition (S,T(S)), and view this as an undirected graph, say H. Assume w.l.o.g. that this 
graph is connected (otherwise we prove the lemma for each connected component of H). Pick 
a spanning tree in H. 

Pick any good j G S, and find a path in the spanning tree from j to each good j' G S. If j 
reaches j' with a path of length 21, then py = apj/b where a and b are products of / utility 
parameters (tijfc's) each. Since alternate edges of this path contribute to a and b, we can 
partition the Ui^s of edges in the spanning tree into two sets, T\ and Ti, such that a uses Ui^s 
from T\ and b uses those from T2. 

Next, consider cti = Ci/ r y i , for i G r(S'). Now, r y i = Uij>/pj', where (i,f) is any edge in the 
network. Find the path in the spanning tree from i to j and use the first edge on this path for 
computing r y i (it is easy to see that all these edges come from set T2), and substitute py using 
the expression stated above, i.e., py = apj/b. 

Since S is a tight set, 
i'es ieT(s) 
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In this equation, substitute for py and q« using the expressions constructed above to get an 
equation with one variable, i.e., pj. Now, it is easy to see that the denominator of pj is < A. 
□ 

Lemma 27 In Stage II, consider two phases P and P', not necessarily consecutive, such that 
good j lies in the newly tight sets at the end of P as well as P' . Then the increase in the price 
of j, going from P to P' , is at least 1/A 2 . 

Proof : Let the prices of j at the end of P and P' be p/q and r/s, respectively. Clearly, 
r/s > p/q. By Lemma 26, q < A and r < A. Therefore the increase in price of j, 

r p J- 
s q ~ A 2 ' 

□ 

Lemma 28 In Stage II, a phase consists of at most g iterations. 

Proof : After each iteration, other than the last one, at least one good will move from G — J 
to J. □ 

The structure of the rest of the argument is quite similar to that of Stage I. Once again, the 
central fact established is that $ drops by an inverse polynomial factor, of (1 — 1/n 2 ), in a 
phase (Lemma 33). Assume that a given phase consists of k iterations. Let Iq denote set / at 
the start of the phase and let // denote the set I at the end of the Z-th iteration, 1 < I < k. 
Assume that at the start of this phase, maxj € s{|^|} = S = Sq. Let 

Si = min{|#j|}, for 1 < I < k, and 9k = 0. 

As in Stage I, we will account for the drop in $ in two steps in each iteration. First, as prices 
of goods in J are increased, the 0j's of buyers i G I decrease, leading to a reduction in <!?. 
Second, when a new edge (j, i), with j £ (G c — J) and i £ I, is added to the network, the flow 
becomes more balanced, leading to a further drop. As in Stage I, we will account for the first 
drop using the l\ norm and the second drop using the I2 norm. 

We begin by accounting for the second decrease. Just before new edge (J, i) is added, let iV 
be the network and p be the prices of goods. Let N' be the network obtained by adding this 
edge to iV; of course, the prices remain unchanged. Let / and /* be balanced flows in TV and 
N', respectively. Denote by (9*) the surplus vector w.r.t. flow / in N (flow /* in N'). 

Lemma 29 \\9\\ 2 - \\9* \\ 2 > 5 2 , where 5 = t - 9*. 

Proof : Since the Invariant holds and the prices are unchanged, / and /* have the same 
value. Therefore, flow f* — f will consist of circulations. Since / is a balanced flow, all these 
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circulations must use the edge (j, i), because otherwise a circulation not using edge could 
be used for making / more balanced. These circulations will have the effect of decreasing the 
surplus of buyer i £ I, and increasing the surplus of buyers i\ G (B — I), for 1 < I < k. Let 
0* - 6 k = 5 h for 1 < I < k. Then, £f =1 5 t = 5. 

For each buyer ii, 1 < %i < k, there is a path from ii to i in the corresponding circulation and 
hence there is a path from i to i\ in the residual graph w.r.t. flow /*. Since /* is balanced, by 
Property 1, the surplus of buyer i w.r.t. /* is at least as large as that of i/. Therefore, 9* > 9* r 
The inequality ||6>|| 2 — \\6* || 2 > 5 2 now follows from Lemma 30, on substituting a = 6*, and for 
1 < I < k, bi = 0*. □ 

Lemma 30 If a > bi > 0, 1 = 1, 2, . . . , k and 5 = Y$=i $l where 5, 5i > 0, 1 = 1, 2, . . . , k, then 
\\(a + S,h-Si,b 2 -S2,...,b k - 5 k ) || 2 - ||(a, h, b 2 , . . . , b k )\\ 2 > 5 2 . 

Proof : 

k k k 

(a + <5) 2 + ]T {h - Si) 2 - a 2 - ][> 2 > 5 2 + 2a(5 - ^ Si) > 5 2 . 

i=l i=l i=l 

□ 

Let 9° denote the surplus vector at the start of the phase and let 6 l denote the surplus vector 
at the end of iteration I, for 1 < I < k. 

Lemma 31 In the l-th iteration, there is a buyer i £ whose surplus decreases by at least 
(<5,_i - Si), forl<Kk. 

Proof : By the definition of set I' in procedure Update sets (II) and Property 1, there is 
a buyer i € I\-\ which achieves minj g / ; {|0j|} at the end of iteration /. Clearly, the surplus of 
i decreases by at least (<5/_i — b~i) in the Z-th iteration. □ 

Lemma 32 For 1 < I < k, 

II^T-II^II 2 ^-!-^) 2 . 

Proof : We first prove the statement for 1 < I < k. By Lemma 31, there is a buyer i E 
whose surplus decreases by at least (<^_i — <5/) in the /-th iteration. Let us split this decrease 
into two parts, the decrease due to increase in the prices of goods in J and that due to a new 
edge entering the network. Let these be a and b, respectively. Clearly, a + b > Si-\ — Si. 

Let 9' be the surplus vector just before the new edge is added to the network in iteration I, i.e., 
right after all the increase in prices of J has happened. As prices in J increase, the surpluses 
of buyers in / decrease, but those of buyers in B — I remain unchanged. Let c be the surplus 
of buyer i at the beginning of iteration /. Then, 

||^- 1 || 2 -||0'|| 2 > c 2 -( c -a) 2 >a 2 + 2ac. 



30 



By Lemma 29, 

\\e'\\ 2 -\\e l f>b 2 . 

Adding the two we get 

||^-l||2 _ ||0l||2 > a 2 + 2ac + b 2 > ( a + ft)2 > ( (Jj _ l _ (5j )2 j 

where the second last inequality follows from the observation that 6 < c. 

Finally, in the /c-th iteration, there is a buyer i G Ik-\ whose surplus changes from 0j > to 
0. Therefore, 

since 5k-i < 6*« and <5fc = 0. □ 
Lemma 33 In a phase in Stage II, the potential drops by a factor of 




Proof : Now, ||#°|| 2 — ||# fc || 2 can be written as a telescoping sum of k terms, each of which is 
the decrease in the potential in one of the k iterations. Lemma 32 gives a lower bound on each 
of these terms. The total lower bound is minimized when each of the differences (<5/_i — 6i) is 
equal. Now using the fact that Sq = 5 and 5k = 0, we get: 

||0O.|2 _[|0fc||2 > ^! 

k 

Finally, since ||#°|| 2 < nS 2 , and by Lemma 28 k < n, we get: 




□ 

Lemma 34 The execution of Stage II requires at most 

O (n 4 (log n + n log U + log C)j 
max-flow computations. 

Proof : By Lemma 33, the potential drops by a factor of half after 0(n 2 ) phases. At 
the start of the algorithm, the potential is at most n. Once its value drops below 1/A 4 , the 
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algorithm requires at most n more phases to compute equilibrium prices. This follows from 
Lemma 26 and Lemma 27. Therefore the number of phases is 



By Lemma 28 each phase consists of n iterations and each iteration requires n max-flow 



Lemmas 25 and 34 give: 

Theorem 35 Algorithms 1 and 2 solve the decision and search versions of Nash bargaining 
game ADNB using 



max-flow computations. 

10 Postmortem and Was Stage I Really Needed? 

As pointed out in Section 4 in [DPSV08], the primal-dual paradigm operates in a fundamentally 
different way in the setting of a rational convex program than in the setting of an integral linear 
program. In the latter setting, in each iteration, it picks an unsatisfied complementary slackness 
condition and satisfies it. On the other hand, in the former setting, the algorithm starts off 
with a suboptimal solution that can be viewed as relaxing a class of the KKT conditions. It 
then tightens these conditions gradually; when they are all fully tightened, the optimal solution 
has been reached. 

Let us first show that this high level picture applies to Stage II of our algorithm as well. 
Consider the situation right after a balanced flow has been computed at any point in Stage II, 
and consider an arbitrary buyer i. At this point, let fi be the flow sent on edge (i, t) by the 
balanced flow; fi is also the money spent by i in the current allocation. Buyer z's available 
money at this point is mi = 1 + Ci/li- Therefore, 



Let Vi be the total utility derived by i from the current allocation and suppose that Xij > 0. 
Then, 



0(n 2 log (A 4 n)) = 0(n 2 (log ra + nlogU + logC)). 



computations for computing a balanced flow. The lemma follows. 



□ 




/• = m - - p. - l = £l _ p.. 



7i = -^- 



Pj fi (-&) + £' 



This yields 
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Substituting for r y i and rearranging we get 

Pj _ u ij 

{-Pi) Vi-Ci 

To summarize, at any point in Stage II, we have ensured the first two KKT conditions and 
relaxed the last two as follows: 

(1) Vj G G : Pj > 0. 

(2) VjGG: Pj >0^ E i6 B^- = l- 
(3') Vie B, Vj € G : -Sfc > 

(4') Vi G B, Vj G G : Xij > = 

Observe that if prices are not feasible, then for some i, /3j > 0. If so, the relaxed KKT 
conditions (3') and (4') will be meaningless. Recall that throughout Stage II, < — j3 i < 1 
and (— monotonically increases and reaches 1 at termination. Thus, at termination, the 
last two KKT conditions are also ensured. For establishing a bound on the number of phases 
needed in Stage II, it suffices to study the potential function 

Clearly, <£' > at the start of Stage II and increases monotonically to n. However, it turns 
out to be more convenient to study the potential function <I> given in Section 9.2, which clearly 
achieves the same end. 

Stage I determines feasibility of the given market. What if we only wanted to solve the 
promise problem of finding the equilibrium of a given feasible market? Since the prices found 
by Initialization are guaranteed to be small, could we not go directly to Stage II and raise prices 
until equilibrium is attained? The answer is "No" . The reason is that Stage II is guaranteed 
to converge only if it is started with feasible prices. As can be seen from the proof of Lemma 
10, if for some buyer i, /3 f > 0, then raising x will actually increase her surplus. This happens 
because her money will increase at a faster rate than the rate at which flow on edge (i,t) 
increases. 

Thus, Stage I not only determines feasibility but, if the given market is feasible, it also finds 
a suitable initial price vector for Stage II. An interesting aspect of Stage I is that it needs to 
decrease prices of goods, even though it starts with a small price vector. Furthermore, it is 
easy to see that if the given market is feasible and Stage I is started off with any small price 
vector, not necessarily the one found by Initialization, it will terminate with a feasible price 
vector. Hence we get the following interesting fact: 

Lemma 36 Let A4 be a feasible flexible budget market and let p be small, positive prices for 
it. Then, there exist positive prices q such that p weakly dominates q and prices q are feasible. 
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Finally we address the question of whether a separate procedure was needed for testing feasi- 
bility, especially in light of the fact that such a procedure is not needed when the primal-dual 
paradigm is used for solving a non-total integral linear program. Let us illustrate the latter by 
comparing, at a high level, the algorithms for the problems of maximum weight matching and 
maximum weight perfect matching in bipartite graphs; for full details, see [CCPS98]. 

The only difference in the LP's of these two problems is that whereas the former demands at 
most 1 matched edge incident at each vertex and the latter demands exactly 1 edge. As a 
result, in the dual, the vertex variables are constrained to be non- negative in the former and 
unconstrained in the latter. This gives rise to an additional complementary slackness condition 
in the former, i.e., any vertex with a positive dual must be matched. 

The algorithm for the former problem attempts to "repair" this complementary slackness 
condition one vertex at a time. It attempts to find an augmenting path from a violating 
vertex, say v, by growing a "Hungarian tree" rooted at v. If another unmatched vertex enters 
the tree, an augmenting path can be found and v is matched off. On the other hand, if the 
tree becomes maximal without encountering another unmatched vertex, then the algorithm 
is able to drive the dual of some vertex in the tree down to zero. If this vertex is v, the 
complementary slackness condition at v has been repaired. If it is some other vertex, say u, 
then there is an alternating path between u and v, which enables the algorithm to unmatch 
u and instead match off v. This repairs the complementary slackness condition at v without 
creating a violation at u. 

The algorithm for the latter problem is very similar - it attempts to iteratively match off 
unmatched vertices. It tries to find an augmenting path from an unmatched vertex, say v, by 
growing a "Hungarian tree" rooted at v. If the tree becomes maximal without encountering 
another unmatched vertex, then the number of "outer" vertices in the tree exceeds the number 
of "inner" vertices by 1. Now, decreasing the duals at outer vertices by A and increasing them 
at inner vertices by A yields a feasible dual. By letting A — > oo, we get that the dual LP is 
unbounded. This gives a proof of infeasibility of the primal LP. 

Consider the two problems ADNB and Fisher's linear case with the money of each agent being 
unit. Observe that the latter problem is a special case of the former when the disagreement 
utilities of all agents are zero, and the latter is total whereas the former is non-total. Why was 
the algorithm for the former so much more elaborate than that for the latter (i.e., the DPSV 
algorithm), especially in view of the fact that the algorithm for maximum weight perfect 
matching is not more involved than that for maximum weight matching? 

In the case of perfect matching, feasibility was ensured one vertex at a time, and when it 
could not be ensured, we got a proof of infeasibility right away. In the case of rational convex 
programs, the KKT conditions were enforced gradually and these conditions were fully satisfied 
only right at the end. So, why can't we simply run Stage II right to the end and then obtain 
a solution or find out that the given instance was infeasible? The reason is that to guarantee 
termination of Stage II, we need to start it off with a feasible price vector, i.e., we need to 
determine feasibility before starting with Stage II. Hence Stage I appears to be essential. 
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11 /i-norm Does Not Suffice 



We give a family of examples showing that the DPSV algorithm, for Fisher's linear case, may 
end up making only inverse exponential progress in a phase if the potential function used is 
the l\ norm of the surplus vector. 

We will define the example in terms of 2 parameters, 5 and H, which will be set at the end. 
Assume B = {&o, &i, • • • , b n -i, b n } and G = {go, g\, . . . , g n -i, 9n}- At the start of the phase, 
the only edges present in the network are {gi, for < i < n. The money of the buyers are 
as follows: 

5 6 

mo = 1 + 5, and for 1 < i < n — 1, mi = — -, and m n = H H — . 

2 l n 

The prices of goods are as follows: 

8 6 
Po = 1 , and for 1 < i < n — 1 , pi = — , and p n = H H . 

2* n 

Hence, at the start of the phase, the surplus of bo is 5, and that of the rest of the buyers is 0. 

We will set 5 = 1 and H to be a large number, say n. The phase starts with / = {60} and 
J = {go}- Assume that at the end of iteration i, edge enters the network, and as a 

result, b{ enters / and gi enters J, for 1 < i < n. The increment in price in each iteration is 
very small - this is easily arranged by choosing the right utilities u^'s. 

To keep the description clean, let us assume the increments in price are all zero; the numbers 
can be easily modified by inverse exponential amounts to yield the desired outcome, even if 
the prices need to increase in each iteration. If so, at the end of all this, the surplus of 6« is 

for < i < n — 1, 



2 i+i' 

and that of b n is 2 ^-i ■ 

Finally, in iteration n + 1, a very slight increase in x leads to set {g n } going tight. Observe 
that the reason for choosing H to be a large number is the ensure that this slight increase in 
x will not make a larger set go tight. Observe that r({g n } = {b n ,b n+ \}. Now, the increase in 
the price of g n needed for this is Since H is a large number and the increase in x is very 

small, the total increase in the prices of other goods is at most a constant factor more. 

In summary, the total increase in the norm of p in this phase is an inverse exponential factor. 



12 Discussion 

Our paper provides two new pieces of evidence to show that the notion of balanced flow is basic 
to the problems tackled in [DPSV08] and the current paper. The first is the result of Section 
11. Second, observe that the notion of a feasible price vector was defined using balanced flows. 
In view of Lemma 6 and the remarks made after it, and of Lemma 36, characterizing the sets of 
small and feasible price vectors for a feasible flexible budget market is an interesting question. 
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As in the case of the EG-program, the optimal solutions of convex program (2) resemble those 
of a linear program rather than a nonlinear program. So, we repeat a question raised in [Vaz07] 
namely, can the solution to ADNB be captured via a linear program? We believe the answers 
to these questions are "no" and that establishing this in a suitable formal framework will 
provide new insights into the boundary between linear and nonlinear programs. 

The most prominent problem for which a rational convex program is known but a combinatorial 
algorithm is not known is the linear case of the Arrow-Debreu market model, i.e., the convex 
program of Jain [Jai07]. We pose the following easier question: Give a polynomial time 
algorithm for this problem that uses only an LP solver. A much more general question along 
these lines is stated in the Introduction. 

We list some more rational convex programs and leave the problem of finding combinatorial 
algorithms for them. First, generalize ADNB to additively separable, piecewise-linear, con- 
cave utilities. This problem has a rational convex program, and the question of finding a 
combinatorial algorithm for it becomes even more significant in view of recent results showing 
PPAD-completeness of computing an equilibrium in the Arrow-Debreu model with these utility 
functions [CDDT09, CT09, VY10]. 

In an interesting paper, Kalai [Kal77] relaxed Nash's axiom of symmetry and derived the solu- 
tion concept of nonsymmetric bargaining games. The convex program capturing the solution 
to the nonsymmetric extension of ADNB is also rational; moreover, this convex program 
generalizes the Eisenberg-Gale program and hence captures Fisher's linear case as well. De- 
spite substantial effort, this problem has not yielded to a combinatorial algorithm. Once 
it is obtained, one could consider the common generalization of the last two problems, i.e., 
nonsymmetric ADNB with additively separable, piecewise-linear, concave utilities. 

On restricting ADNB (nonsymmetric ADNB) to zero disagreement utilities we get the prob- 
lems of computing equilibria for linear Fisher markets with unit (arbitrary) money among 
buyers. Of course, both these problems are total. It turns out that a combinatorial algorithm 
for the unit money case is no easier than that for the arbitrary money case. In view of this, the 
difficulty of obtaining a combinatorial algorithm for nonsymmetric ADNB comes as a surprise 
and may be substantiating the observation that in the setting of rational convex programs, 
non-total problems behave quite differently from total problems. 

All the rational convex programs mentioned above involve the log function in the objective, 
together with linear constraints. Another class of rational convex programs is obtained by 
having a quadratic objective function and linear constraints. It will be very interesting to 
obtain combinatorial polynomial time algorithms for such rational convex programs as well. 

The reader can see that the algorithm for ADNB exploits a surprisingly rich and clean struc- 
ture which is, in some ways, reminiscent of the majestic structure of matching. In our experi- 
ence, such structure does not occur in isolation and we believe that what we see so far is the 
tip of an iceberg. This leads to the question: what does the rest of the iceberg look like? 

One possibility is to seek combinatorial approximation algorithms for solving specific classes 
of nonlinear convex programs. In this respect, important hints may be obtained from the 
way the primal-dual paradigm was extended from solving linear programs exactly to obtaining 
near-optimal solutions to linear programs within the area of approximation algorithms. The 
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mechanism involved in all of the latter algorithms was that of relaxing complementary slack- 
ness conditions, which was first formalized in [WGMV95]. We also note that in the setting of 
approximation algorithms, the primal-dual paradigm has been successful primarily for mini- 
mization problems. So, our more precise question is, "Is there a natural way of relaxing the 
KKT conditions to obtain primal-dual (combinatorial) algorithms for near-optimally solving 
interesting classes of (perhaps minimization) nonlinear convex programs?" 
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A Solution to ADNB in the Limit 

Assume that we are given an instance of game ADNB that is feasible and let M be the flexible 
budget market obtained from it. In this section, we present an algorithm that converges to 
the equilibrium of M in the limit. 

Algorithm 3 will use the DPSV algorithm as a subroutine [DPSV08]. When this subroutine is 
called, we assume that the money of each agent is fixed and is specified in the vector m. 

Let N'(p) denote the network for the case that the money of agents is fixed and specified by 
vector m; this network differs from network N(p) only in that the capacities of edges going 
from buyers to t are specified by m, rather than being defined as a function of the prices. 



Algorithm 3 (Solution to ADNB in the Limit) 

1. Initialization: Vi G B : mj <— 1. 

2. Compute equilibrium prices, p, for market (it, m) using the DPSV algorithm. 

3. For each i E B, compute r ) i w.r.t. prices p, and set m\ «— 1 + 

4. If m! = m then output equilibrium allocations and HALT. 
Else, update m to m' and go to Step 2. 
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Let p* and m* be the equilibrium prices and moneys for the flexible budget market M , and let 
p( fe ) and m( fc ) denote the prices and moneys computed by the algorithm in the k-th iteration, 
k > 1. 

Lemma 37 p^ and are monotone increasing and are weakly dominated by p* and m* , 
respectively. 

Proof : We will use the following 2 facts. First, the DPSV algorithm maintains the follow- 
ing invariant throughout: 

Invariant: W.r.t. current prices, p, (s, B U G U t) is a min-cut in network N'(p). 
Second, if p are equilibrium prices for money m and if m' is at least as large as m in each 
component, then the equilibrium prices for money m' cannot be smaller than p in any com- 
ponent. 

Consider the following induction hypothesis: 

• the algorithm given above maintains the Invariant throughout. 

• p( fc ) is monotone increasing (hence, for each agent i, 7 i is monotonically decreasing). 

• m,( fe ) is monotone increasing. 

It is easy to carry out this induction simultaneously for all 3 assertions. 

Using the first assertion and Lemma 6, p^> is weakly dominated bounded by p* . Now, using the 
formula for money in flexible budget markets, it is easy to see that is weakly dominated 
by m*. □ 

Theorem 38 Algorithm 3 converges to the equilibrium prices and moneys of market M. in 
the limit. 

Proof : We will use the following fact: for the linear case of Fisher's model, the analog 
of Lemma 4 holds, i.e., if p are equilibrium prices for money m, then in network N'(p), 
(s,BL)GUy) and (sUBLiG,t) must both be min-cuts (for a proof, see Lemma 5.2 in [Vaz07]). 

Since p^ and are monotone increasing and bounded, they must converge. Let p^ and 
be their limit points. W.r.t. these prices and moneys, it must be the case that for each 
i G B, rrii = 1 + Ci/j i and (s, B U G U t) and (sUBUG,t) must both be min-cuts in the 
corresponding network (by the fact stated above). Using lemma 4 we get that p^ -* and 
are equilibrium prices and moneys for market M.. □ 

Finally, by Theorem 3 we get: 

Corollary 39 Algorithm 3 converges to the Nash bargaining solution for ADNB. 
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